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INTRODUCTION 

Nucleophilic  substitution  has  long  been  recognized  as  one  of  the  most  important  reactions 
in  synthetic  organic  chemistry.  Because  of  their  predictable  regio-  and  stereochemistry, 
reactions  proceeding  through  the  SN2  mechanism  have  received  considerable  attention.  The  ease 
of  displacement  of  the  leaving  group  plays  a  significant  role  in  determining  the  reactivity  of  a 
given  substrate;  however,  even  compounds  that  contain  extremely  poor  leaving  groups  such  as 
ethers,  amines  or  thioethers  are  reactive  and  synthetically  useful  if  these  functional  groups  are 
structured  in  such  a  manner  that  the  heteroatom  is  incorporated  into  a  three-  or  four-membered 
ring.  Relief  of  ring  strain  in  the  transition  state  is  typically  cited  as  the  sole  source  of  the 
increased  reactivity  of  heterocycles  with  nucleophiles  relative  to  acyclic  analogs. 

Over  a  decade  ago,  however,  it  was  noted  that  ring  strain  alone  is  insufficient  to  account 
entirely  for  rate  increases  in  SN2  reactions  relative  to  suitably  chosen  acyclic  model  compounds. 
Lillocci  found  the  rate  of  the  cleavage  reaction  of  an  aziridinium  triflate  with  acetonitrile  in  the 
presence  of  N-ethyldiisopropylamine  to  be  at  least  103  times  faster  than  that  of  the  corresponding 
azetidinium  salt  in  spite  of  similar  ring  strain  energies.1 2  This  phenomenon  has  also  been 
observed  in  nucleophilic  cleavage  of  these  ions  by  sodium  methoxide.3  In  a  related  reaction, 
Stirling  observed  that  the  rate  of  ring  cleavage  of  cyclopropanols  under  basic  conditions  was 
considerably  faster  than  that  of  cyclobutanols.4  Recently,  Hoz  recognized  that  an  “additional 
factor,”  as  yet  uncharacterized,  must  be  included  to  explain  the  high  computed  reactivity  of 
anionic  nucleophiles  with  3-membered  relative  to  4-membered  heterocycles.5  Although  Houk 
and  coworkers6  have  recently  proposed  that  orbital  interactions  through-bonds  are  responsible  for 
the  large  rate  enhancement  in  cleavage  of  three-  with  respect  to  four-membered  rings,  the  present 
investigation  was  designed  to  discover  an  alternative  explanation  of  this  rate  increase  in  terms  of 
disjoined  effects,  present  for  three-  but  not  four-membered  rings  after  strain  relief  in  the  transition 
state  has  been  taken  into  account. 

We  have  been  interested7'8  in  the  studying  relative  rates  of  the  reactions  of  small  ring 
heterocycles  with  amines.  These  Menschutkin-type  reactions  are  examples  of  the  general  type  of 
Sn2  reactions  in  which  an  uncharged  nucleophile  reacts  with  an  uncharged  substrate.  Charge 
separation  occurs  along  the  reaction  coordinate  as  the  nucleophile  assumes  a  partial  positive 
charge  and  the  heteroatom  a  partial  negative  charge.  Due  to  extreme  endothermic ity  occasioned 
by  separation  of  developing  opposite  charges  in  the  transition  state  in  the  absence  of  stabilizing 
intereactions  with  solvent,  for  example,  attempts  to  gain  a  fundamental  understanding  of  this 
reaction  in  the  gas  phase  have  been  limited  to  one  experimental  study.9  This  reaction  type  can  be 
contrasted  with  the  extensively  studied10  gas  phase  “anionic”  SN2  reaction,  exemplified  by  the 
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attack  of  bromide  ion  on  methyl  chloride.  Because  Menschutkin-type  reactions  are  of 
considerable  synthetic  importance  in  solution,11  a  fundamental  understanding  of  this  chemistry  in 
the  absence  of  potentially  complex  solvent  effect  is  desirable.12  Computational  studies 
combining  an  appropriate  method  and  basis  set  of  reasonable  size  are  a  cost-effective  approach  to 
elucidating  the  mechanistic  parameters.13 

Oxirane  and  oxetane  have  nearly  the  same  ring  strain  energies:  27.3  and  25.5  kcal/mol, 
respectively,  while  for  thiirane  and  thietane  these  values  are  19.8  and  19.6  kcal/mol, 
respectively.14  These  pairs  of  compounds  offer  an  opportunity  to  address  the  question  of  whether 
ring  strain  relief  is  largely  responsible  for  their  rates  of  reactions  with  nucleophiles.  It  would  be 
expected  that  if  ring  strain  release  in  the  transition  state  occurs  to  equal  degrees  in  each  pair  of 
compounds,  the  rates  of  reaction  would  be  very  similar.  In  this  investigation,  we  undertook  a 
computational  investigation  of  the  rates  of  reaction  of  each  of  these  compounds  with  ammonia  in 
the  gas  phase.  In  the  absence  of  solvent  effects,  the  results  may  be  interpreted  in  terms  of  the 
intrinsic  properties  of  the  reaction.  In  order  to  assess  the  effect  of  ring  strain,  the  rate  of  reactions 
with  acyclic  model  compounds  were  also  investigated. 


COMPUTATIONAL  METHODS 

The  computations  were  performed  in  the  gas  phase  using  Gaussian  94. 15  Clearly  a 
method  that  incorporates  electron  correlation  is  required  for  these  studies  of  nucleophilic 
displacement  reactions.  The  6-31+G(d)  basis  set  which  includes  diffuse  functions  to 
accommodate  the  lone  pairs  and  anionic  species  found  in  these  reactions  was  chosen.  Two 
methods,  B3LYP  and  MP2(Full),  were  evaluated  in  terms  of  their  ability  to  reproduce  the 
reported  strain  energies  of  the  heterocycles  studied  using  the  method  of  Dudev  and  Lim. 16  Zero 
point  correction  factors  of  0.9804  and  0.9646  were  used,  respectively.17  As  can  be  gleaned  from 
Table  1,  the  MP2(Full)  method  agreed  more  closely  with  the  experimental  values,  and  was 
selected  for  the  subsequent  rate  calculations. 


Table  1.  Strain  Energies  and  Verification  of  Computational  Method  (kcal/mol) 


Experimental14 

MP2(Full)/6- 

31+G(d)//MP2(Full)/6- 

31+G(d) 

B3LYP/6-3 1+G(d)// 
B3LYP/6-3 1+G(d) 

Oxirane 

27.3 

26.1 

31.2 

Thiirane 

19.8 

19.3 

20.8 

Oxetane 

25.5 

26.3 

28.8 

Thietane 

19.6 

23.1 

22.1 

Transition  states  and  ground  states  were  identified  by  the  presence  of  one  and  zero 
imaginary  frequencies,  respectively.  Transition  states  were  further  corroborated  by  animation  of 
the  imaginary  frequency  using  GaussView  2.03 18  and  by  means  of  intrinsic  reaction  coordinate 
(IRC)  calculations  to  illustrate  that  the  calculated  structure  was  found  on  the  potential  energy 
surface  connecting  the  reactants  and  the  products. 

The  strain  energy  of  a  cyclic  molecule  was  measured  relative  to  a  suitable  acyclic  analog, 
methoxyethane  or  thiomethylethane.  The  most  stable  conformation,  initially  obtained  by  means 
of  an  AMI  semi-empirical  conformational  search,  was  used  for  computations  for  the  acyclic 
compounds.  One  of  the  difficulties  inherent  in  calculating  relative  rates  of  cyclic  with  respect  to 
acyclic  compounds  is  that  there  are  no  perfect  acyclic  models  for  cyclic  compounds.  Dimethyl 
ether  and  dimethyl  thioether  were  rejected  as  models  even  though  they  possess  the  correct 
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number  of  carbon  atoms  as  oxirane  and  thiirane  because  they  lack  secondary  carbons.  The 
reactions  of  ammonia  at  the  secondary  carbons  of  methoxyethane  and  thiomethylethane  were 
deemed  the  best  choices  to  model  reactivity  for  the  oxygen  and  sulfur  heterocycles,  respectively. 

Calculations  using  atoms  in  molecules  theory  (  AIM)  were  performed  by  means  of  the 
AIM  2000  program19 


DISCUSSION  AND  RESULTS 

Relative  rates  of  reaction  with  ammonia  for  each  series  of  compounds  are  provided  in 
Table  2.  At  this  level  of  theory,  IRC  calculations  identified  dipolar  complexes  for  oxirane, 
thiirane,  oxetane  and  thietane  that  were  3.2,  2.6,  0.3  and  2.2  kcal/mole,  respectively,  more  stable 
than  the  separated  reactants;  no  species  other  than  the  separated  reactants  could  be  identified  for 
the  remaining  (acyclic)  substrates 

(Table  S-l.)  These  small  differences  in  stability  may  be  due  to  the  fact  that  the  heterocycles  have 
their  alkyl  groups  pinned  back  so  as  to  avoid  steric  interaction  with  ammonia  when  it  approaches 
close  enough  to  produce  an  effective  dipole-dipole  interaction.  Not  unexpectedly,  both  three- 
and  four-membered  heterocycles  were  very  much  more  reactive  than  their  corresponding  model 
compounds.  Moreover,  the  3-membered  rings  wee  more  than  106  times  more  reactive  than  the  4- 
membered  heterocycles. 

In  order  to  address  disjoined  effects,  the  amount  of  strain  energy  remaining  at  the 
transition  state  must  be  determined.  The  total  strain  of  a  molecule  in  the  gas  phase,  V,  is  given 
by20 


V  V stretch  ^  V bend  T  V torsion  T  V  VDW 

where  the  sums  of  bond  stretching  and  bending  energies  and  the  sums  of  torsional  and  nonbonded 
interactions  are  included.  Stretching  energies  and  nonbonded  interactions  are  unimportant  for 
these  transition  state  structures.  For  purposes  of  this  discussion,  it  will  be  assumed  that  the 
bending  energies  and  torsional  interactions  involving  the  atoms  of  the  reaction  center  are  small 
and  approximately  equal  for  the  members  of  each  series.  Inspection  of  the  transition  state 
geometries  establishes  that  this  is  a  reasonable  supposition.  Using  the  usual  harmonic  potential 
based  on  Hooke’s  Law  functions  with  0.025  kcal  mol'1  deg 2  as  the  constant  for  angle  bending, 21 
and  using  3.5  kcal  mol'1  and  3.3  kcal  mol'1,  the  torsional  barriers  for  ethanol22  and  ethanthiol,23 
respectively,  for  the  torsional  barriers  in  the  corresponding  transition  states,  one  can  estimate  of 
the  sum  of  the  angle  deformation  and  torsional  energies  of  the  transition  states  for  oxirane, 
oxetane,  thiirane  and  thietane  to  be  3.3,  4.6,  2.7  and  1.5  kcal  mol’1,  respectively.  If  these 
estimates  are  subtracted  from  the  calculated  strain  energies  of  the  heterocycles,  the  strain  energy 
released  in  the  transition  state  for  each  substrate  is  obtained.  Addition  of  these  values  to  the 
enthalpies  of  activation  (Table  2)  for  each  reaction  permits  one  to  calculate  G*  =  69.3,  77.6, 

60.7  and  74.7  kcal  mol’1,  respectively  for  oxirane,  oxetane,  thiirane  and  thietane  if  each 
contribution  from  strain  relief  were  removed.  The  calculated  values  for  the  4-membered  rings 
agree  quite  well  with  the  expected  values  based  solely  on  relief  of  ring  strain  relative  to  the  strain- 
free  acyclic  model  compounds  (  G±=  78.2  and  74.4  kcal  mol'1  for  MeOEt  and  MeSEt, 
respectively.)  The  3-membered  rings  are  quite  exceptional.  Including  the  correction  for  strain 
energy,  oxirane  is  found  to  react  with  ammonia  at  a  rate  3.27  x  106  times  faster  than  methyl  ethyl 
ether.  Thiirane  reacts  at  a  rate  of  1.08  x  1010  faster  than  its  acyclic  model  compound. 

The  structure  of  the  transition  state  is  critical  to  understanding  the  relative  rates  of  these  reactions. 
While  the  geometry  of  the  transition  states  is  useful,  it  fails  to  address  directly  the  distribution  of 


3 


electrons  in  these  species.  Bader’s24  theory  of  atoms  in  molecules  (AIM)  is  ideally  suited  to 
studying  the  electronic  distribution  of  the  transition  state.  This  approach  has  demonstrated  that 
the  negative  of  the  Laplacian  of  the  electron  density  coincides  with  electronic  charge 
concentration  (bonding)  in  treatments  of  localized  electron  density  by  the  Lewis  and  VSEPR 

—  2 

models.  The  results  of  an  AIM  calculation  of  -V  P  for  a  typical  reaction,  attack  of  ammonia  on 
thiirane,  are  provided  in  Figure  1.  This  relief  graph  clearly  demonstrates  the  localization  of 
electron  density  between  the  bonded  atoms,  i.e.,  the  carbons  and  between  sulfur  and  the  carbon 
that  is  not  under  attack  by  ammonia.  It  is  noteworthy  that  electron  density  is  absent  in  the  region 
between  the  reaction  center  and  sulfur,  while  capture  of  electron  density  from  nitrogen  is  only 
modest.  Some  years  ago,  doubt  was  cast  on  the  classical  representation  of  the  transition  state  in 
Sn2  reactions  that  involves  essentially  equivalent  degrees  of  bond  making  and  bond  breaking. 
Bader,  Duke  and  Messer25  calculated  that  the  SN2  reactions  of  fluoride  or  cyanide  ion  with 
fluoromethane  proceed  through  a  partially  positively  charged  carbon  fragment.  The  present 
results  indicate  that  Menschutkin-type  reactions  also  depart  from  classical  Sn2  mechanistic 
formalism  with  bond  breaking  significantly  advanced  relative  to  bond  making.  The  virtually 
complete  cleavage  of  the  carbon-heteroatom  bond  in  the  transition  state  prior  to  formation  of  the 
carbon  nitrogen  bond  requires  that  charge  be  distributed  over  the  carbon  chain.  Molecules  that 
are  more  efficient  lowering  the  energy  of  this  species  will  react  faster. 


Figure  1.  Faplacian  of  the  Electron  Density  of  the  Transition  State  for  Thiirane  +  NH3 
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Table  2.  Relative  Rates  of  the  Reactions  of  Heteocycles  and  Acyclic  Analogs  with  NH 
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What  effect  or  effects  are  producing  such  dramatic  rate  increases?  Since  separation  of 
opposite  charges  occurs  along  the  reaction  coordinate,  it  appeared  reasonable  to  calculate  atomic 
charges  in  the  transition  states.  Small  changes  in  charge  can  produce  profound  changes  in 
energy.26  If,  for  example,  a  charge  of  +0.237,  separated  at  a  distance  of  2.08  A  from  a  charge  of  - 
1 . 1 3,  is  increased  to  +0.250,  the  energy  of  the  system  along  the  reaction  coordinate  required  to 
maintain  separation  increases  by  1 1.3  kcal  mol"1. 

Interpretation  of  the  results  of  atomic  charge  calculations,  however,  must  be  approached 
with  caution.  Atomic  charges  are  not  physical  observables;  they  are  dependent  upon  the  method 
by  which  they  are  derived.  The  calculated  charges  are  not  necessarily  centered  on  atomic 
nuclei.27  An  excellent  discussion  that  compares  various  methods  for  calculating  atomic  charges 
is  available.28 

The  approach  chosen  here  was  to  calculate  the  attractive  energy  using  the  atomic  charges 
obtained  from  four  methods:  Mullikan,  AIM,  CHELPG,  and  NBO.  The  results  are  provided  in 
Table  3  with  the  charges  of  the  hydrogens  summed  into  the  heavy  atoms  to  which  they  are 
attached.  Using  a  classical  electrostatic  model,  a  dielectric  constant  of  1.0, 29  and  neglecting  the 
charges  on  all  atoms  other  than  the  reaction  center  and  leaving  group,30  the  attractive  energy 
difference  between  oxetane  and  oxirane  is  20  kcal  mol'1,  and  that  for  thietane  and  thiirane  is  12 
kcal  mol'1,  using  AIM  atomic  charges.31  The  higher  the  attractive  energy,  the  greater  the  energy 
required  to  reach  the  transition  state.  These  estimated  energies  are  more  than  sufficient  to 
account  for  rate  accelerations  of  oxirane  and  thiirane  above  and  beyond  that  due  to  relief  of  ring 
strain.  The  same  conclusion  is  reached  using  the  values  obtained  using  the  remaining  atomic 
charge  methods  with  the  sole  exception  the  NBO  method  for  the  sulfur  heteocycles. 

A  reasonable  explanation  for  the  charge  distributions  observed  in  these  transition  states 
for  the  ethers  was  suggested  by  Wiberg32  in  his  discussion  of  calculated  results  for  alkoxide  ions 
at  the  MP3/6-3 1 1++G**  level.  The  negatively  charged  oxygen  was  postulated  to  repel  electron 
density  of  the  back  lobe  of  the  adjacent  atom,  resulting  in  the  adjacent  atom  becoming  relatively 
positive  and  the  atoms  attached  to  it  relatively  negative  with  respect  to  the  corresponding  alcohol. 
To  ascertain  if  these  results  still  obtain  at  our  lower  level  of  theory,  atomic  charge  calculations 
were  repeated  at  the  MP2(Full)/  6-3 1+G(d)  level.  The  results  of  these  calculations  are  presented 
in  Tables  S-2  and  S-3,  proving  that  the  same  interpretation  may  also  be  applied  to  the  present 
calculations. 

A  rationalization  similar  to  that  used  above  for  the  transition  states  for  the  ethers  cannot 
be  applied  to  the  transition  states  for  the  reactions  of  the  thioethers  with  ammonia.  If  the 
distribution  of  charges  in  ethanethiolate  ion  are  calculated  at  the  MP2Full)/6-31+G(d)  level,  it  is 
found  that  negative  charges  reside  on  both  the  geminal  and  vicinal  carbons,  differing  from  the 
results  for  ethoxide  ion  where  large  positive  charge  is  found  on  the  geminal  carbon.  Using 
perturbation  molecular  orbital  theory,  the  destabilizing  ethyl  7t  -  sulfur  3p  interaction  was  found 
to  be  more  important  that  the  stabilizing  ethyl  n*  -  sulfur  3p  interaction.  This  interaction  results 
in  the  repulsion  of  electrons  by  the  negatively  charged  sulfur  and  an  accumulation  of  negative 
charge  on  the  neighboring  carbons. 33  Partial  negative  charge  is  found  at  the  positions  next  to 
sulfur  in  the  transition  states  using  every  method  (  with  the  exception  of  CHELPG),  mirroring  the 
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Table  3.  Atomic  Charges3  Calculated  at  the  MP2(Full)/6-3 1+G(d)  Level  by  Various  Methods 

Mullikan  AIM  CHELPG  NBO  Mullikan  AIM  CHELPG  NBO 
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'Hydrogens  summed  into  heavy  atoms  to  which  they  are  attached 


results  for  ethanethiolate  ion.  It  is  most  interesting  that  the  charge  at  the  position  vicinal  to  sulfur 
for  the  thiirane  transition  state,  which  is  also  the  reaction  center,  is  the  least  positive  of  all 
reaction  centers.  When  this  position  is  compared  to  that  of  the  transition  state  of  thietane  whose 
alkyl  chain  is  one  carbon  longer,  it  is  found  that  there  is  slight  negative  charge  at  this  position. 
This  result  parallels  that  of  the  thiolate  ion,  and  suggests  that  negative  charge  is  also  released  to 
this  position  for  the  thiirane  transition  state.  The  net  result  is  that  the  reaction  center  for  the 
thiirane  transition  state  is  rendered  less  positive  than  the  other  members  of  this  series.  In  the  case 
of  thietane,  the  reaction  center  is  one  carbon  further  removed  along  the  chain  and  unable  to 
significantly  benefit  from  this  electron  release.  For  thiomethylethane,  the  chain  necessary  for  this 
electron  repulsion  mechanism  to  reduce  the  partial  positive  charge  of  the  reaction  center  has  been 
ruptured,  rendering  stabilization  of  its  reaction  center  by  this  process  impossible.  A  similar 
explanation  can  be  readily  applied  to  the  ethers. 

Within  the  assumptions  of  transition  state  theory,34  one  can  calculate  that  thiirane, 
thietane  and  thiomethylethane  react  57.4,  608  and  1 12  times  faster  than  oxirane,  oxetane  and 
methoxyethane.  The  higher  reactivity  of  the  thioethers  is  consistent  with  the  significantly  lower 
atomic  charges  of  the  thioether  transition  states  leading  to  a  decrease  in  the  energy  required  to 
reach  the  transition  state(  Table  3.)  The  simplest  explanation  of  the  larger  charges  for  the  ethers 
is  that  there  is  a  large  difference  in  electronegativity  for  carbon  and  oxygen,  while  the  that  of 
carbon  and  sulfur  is  close  to  zero.35 

Supporting  this  explanation  of  the  accelerated  rate  of  reaction  of  the  three-  with  respect  to 
four-membered  heterocycles  in  nucleophilic  reactions  with  ammonia  is  the  observation  that  the 
calculated  carbon-carbon  bonds  connecting  the  reaction  center  to  the  chain  in  the  transition  states 
for  thiirane  and  oxirane,  1.476  and  1.492  A,  are  significantly  shorter  than  the  corresponding 
distances  for  the  other  transition  states  (Table  4).  These  stronger  bonds  are  occasioned  by  the 
higher  electron  density  found  in  these  bonds  due  to  electron  repulsion  by  sulfur  and  oxygen  in  the 
manner  described  above  so  as  to  decrease  charge  at  the  reaction  center. 

It  should  be  borne  in  mind  that  it  is  not  possible,  given  the  assumptions  in  assigning 
atomic  charges  (vide  infra),  to  quantitate  the  magnitude  of  this  electrostatic  effect.  Other  factors 
may  contribute  to  the  disjoined  effect;  however,  it  is  likely  that  the  major,  if  not  sole  contributor, 
is  the  one  having  an  electrostatic  origin. 

The  results  suggest  that  a  significant  increase  in  reaction  rate  could  be  observed  for  any 
compound  that  is  able  to  attain  a  transition  state  similar  to  those  for  the 
3-membered  heterocycles  in  this  investigation.  In  general,  any  substrate  that  contains  a 
negatively  charged  substituent  at  the  vicinal  position  in  the  transition  state  should  react  at  an 
increased  rate.  In  a  trivial  extension,  one  would  expect  nucleophilic  cleavage  of  aziridines  to 
proceed  considerably  faster  than  azetidines.  Moreover,  it  should  be  possible  to  observe  enhanced 
rates  for  nucleophilic  substitution  reactions  of  acyclic  compounds. 

Before  addressing  this  question  it  seemed  appropriate  to  address  the  question  of  whether 
the  Menshutkin  reaction  for  acyclic  substrates  proceeded  through  a  transition  state  featuring 
significant  rupture  of  the  bond  between  carbon  and  the  leaving  group  as  was  found  for  the  anionic 
reaction  studied  by  Bader.25  This  type  of  transition  state  would  encompass  the  electrostatic 
requirements  described  above.  For  the  reaction  of  CH3C1  and  NH3  studied  by  computationally  by 
Wiberg,13a  the  Laplacian  of  the  electron  density  was  determined  and  shown  to  clearly  indicate  that 
the  transition  state  had  virtually  complete  C-Cl  bond  breaking  with  incomplete  C-N  bond  making. 
Future  research  will  be  directed  toward  study  of  disjoined  effects  with  suitably  chosen  acylic 
model  compounds. 
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Table  4.  Atomic  Charges  for  Ethanethiol :  Cl  C'CTE  SH 


Mullikan 

AIM 

CHELPG 

NBO 

Cl 

-0.0463 

-0.1301 

0.1248 

-0.1110 

c2 

0.0297 

-0.0959 

0.0545 

0.0331 

s 

0.0166 

-0.0607 

-0.1793 

0.0779 

Table  5.  Atomic  Charges  for  Ethanethiolate  Anion  :  2CH31CH2  S 


Mullikan 

AIM 

CHELPG 

NBO 

c, 

-0.0731 

-0.1406 

0.2245 

-0.6174 

c2 

-0.0313 

-0.1118 

-0.1355 

-0.0551 

s 

-0.8956 

-0.7450 

-1.0889 

-0.7321 

CONCLUSIONS 

The  results  presented  above  demonstrate  that  the  accelerated  rate  of  4-membered 
heterocycles  with  ammonia  relative  to  acyclic  model  compounds  was  due  almost  entirely  to 
release  of  ring  strain  in  the  transition  state.  Three-membered  rings  were  calculated  to  react  more 
than  106  times  faster  than  4-membered  rings.  These  acceleration  are  due  to  a  combination  of 
relief  of  ring  strain  and  the  disjoined  effect.  The  best  explanation  of  this  effect  appears  to  be  the 
favorable  topography  for  minimizing  the  attractive  interaction  between  the  partially  negatively 
charged  heteroatom  and  the  reaction  center  that  is  significant  only  in  the  transition  states  for 
three-membered  rings.  Using  different  mechanisms,  electron-rich  oxygen  and  sulfur  can  release 
electron  density  to  the  reaction  center,  minimizing  its  absolute  positive  charge,  and  therefore, 
making  separation  of  the  heteroatom  easier  due  to  decreased  electrostatic  attractive  interactions. 
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